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Dear DUG Members, 

Late, very late indeed, I present 
DNL#96. Due to many reasons, too many 
to explain, I was not able to finish this 
DNL earlier. Please apologize and enjoy 
this issue. 

I'd like to express my gratefulness to so 
many wishes for Christmas and New 
Year. It was impossible to answer to all 
mails, many thanks, we were really 
moved. 

What concerns DNL#96, 1 
invite you to notice the 
announcement of 

ACA2015, held in Greece 
this summer. This is a 
great conference series 
focussing on Computer 
Algebra Applications. It is 
a long tradition to have an 
Educational Session among 
other sessions presenting CA applica- 
tions in various fields. I attended ACA 
several times and it was always exciting 
to learn what can be done with CAS - 
not only in Teaching Mathematics and in 
Pure Mathematics - not to mention the 
high quality lectures in the Educational- 
session. 

Take the chance to combine visiting one 
of the cradles of European Culture 
(Mykene, Sparta, Epidaurus, Olympia, ...) 
and learning about modern mathematics 
applications. 

We have more Nspire contributions in 
this DNL, some of them based on earlier 
DERIVE applications. I received a couple 


of mails complaining: “Oh, how easy was 
this with our good old DERIVE!" 

Thanks to Heinrich Ludwig we can offer 
also a pure DERIVE article. His stitch 
patterns based on complex numbers in- 
vite for further investigations and dis- 
coveries. It would need another LUA- 
program to transfer his patterns on the 
Nspire machine. Please notice his bundle 
of comments and questions at the end 
of his article. Heinrich produced his 
paper some time ago 
assuming that DERIVE's 
MOD- and MODS-function 
don't support complex 
division. It seems to be 
that later versions of 
DERIVE are able to apply 
MOD and MODS as well for 
complex numbers. 

I finish with my best regards until next 
DNL. Btw, we should reach DNL#100 
this year. 

Always Yours 

Josef 



The late Christmas Star was created by 
Roland Schroder for the DUG- 
Community. 



Download all DAX-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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I 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE , 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 


Next issue: March 2015 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encryption, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the Tl, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

A New Approach to Taylor Series, D. Oertel, GER 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 

Rational Hooks, J. Lechner, AUT 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 
Technical Problems solved with Secondary Maths, W. Alvermann, GER 
Pickover’s Mygalomorphs and Spiders, A. Roulier & J. Bohm, SUI/AUT 
Statistics of Shuffling Cards, Charge in a Magnetic Field, H. Ludwig, GER 
What if General Lee Had a Graphing Calculator at Gettysburg? R. Hawkes, USA 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Invitation from Jean- Jacques Dahan 

Enjoy these 3 videos 

[22/01/15 17:07:21] DAHAN Jean-Jacques: http://youtu.be/cszoiCBU13c 
[22/01/15 17:07:49] DAHAN Jean-Jacques: http://youtu.be/OiUHplMRSTk 
[22/01/15 17:08:13] DAHAN Jean-Jacques: http://youtu.be/_WW62xi_TIQ 

Amities 
Jean-Jacques 

Hubert Langlotz, Germany 

Problems with equations 

Dear Josef, 

Working on one of our examination problems 
(Modeling with a pear-formed body) we came 
across a problem with TI-Nspire. 

Partial problem f) reads: 

f{x) = 4 • Vx • e~°' 5x and g(x) = Ja-x + b shall be 
connected in point P so that we have a smooth tran- 
sition (= differentiable). Find a and b. 

The procedure is well known and leads to a nonlinear system of two equations with two unknowns, a 
and b. 



TI-Nspire seems to be unable to solve this system. How easy was this with good old DERIVE: 
#1: f(x) := 4-^/x- EXP(- 0.5-x) 

#2: g(x) := 7( a-x + b) 

#3: SOLVEff (4) = g(4) a f f (4) = g'(4), [a, b]) 

-4 -4 

#4: a = - 4S-e a b = 256-e 


= -0.87915 06666 a b = 4.688803555 



Do you have any idea how to overcome this deficiency? Regards Hubert 
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Hallo Hubert, 

I played a little bit around with your pear. It works - but Nspire needs some support: 

It is very interesting that we get the exact result solving the system for {x,b} - but not for 
{x,a}!! Don’t ask me why? Is x here something like a catalyst? 

I recommend to add a suitable estimation value for one of the unknowns. 

I must agree with you: this does not replace the functionality of DERIVE. 


& 


solve! j^ 4 ) , 


v=C5 and a= 6- (Ji 
solve! j^ 4 ^ 4 ^ . :{x,b} 

solvefj^ 4 )?^, , ,{.x,b}\ 

H/j(4)=*l(4) } 

& solvefM 4 ^ ,^ 4 ^ , ,{a,£>=l} 

lVl(4)= g j(4) [ 

& solve[|^^ v^ 4 ^, , ,{«=l,fc}j 
!V/(4)= g ;(4) I 


b ■ e +144 +12 ■ e and , 24- Jb ■ e +144 +b- e +288 =8 or ,v=CJ and^ 


-4 -4 

x=C6 and a="48- e and £>=256- e 


.v=c7 and «=~0. 879151 and £?=4.6888 


a=-0 . 879151 and fc=4.6888 


a=-0. 879151 and £>=4.6888 


Hubert sent a second example where Tl-Nspire didn’t behave as expected. It is a 7 th order 
equation. The real bad fact is not that Nspire cannot solve the equation but that we receive 
two different solutions. Students will ask: which one is the correct one? 



We have the choice!! 

DNL: I tried with V200, too. nSolve gives the correct solution immediately. In Exact Mode we 
get no solution with solve (that’s clear) but solve approximated gives 1 .8919 - which is a 
third possible solution? 
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I recommend either applying nSolve (because we know that there is no exact solution) or 
working with solve and similar to the problem above add an estimation for the solution. 


^ solve((9-5- _t) 7 -2=0,.t=l 

| 1=1.57918 

solvet(9-5- .r) 7 -2=0,.v=0, 

| *=1.57918 

nSolve((9-5- .r) 7 =2,.r) 

1.57918 

nSolve((9-5-.t) 7 -2=0,.r) 

1.57918 


See again the difference: 


solve((9-5- _r) 7 -2=0,_tj 

r=2. 05738 

solve((9-5- .v) 7 -2=0,.v=l) 

v=1.57918 


And what is always i good idea: check the solution - or at least find an estimation value - in 
the graphic screen: 



Christoph Kuderli from Switzerland sent an answer on my request posed in the revised 
DNL#33 how to produce the Snowflake Curve in a recursive way: 

„I wrote a Koch curve program which does not work in the Graphs-coordinate system but in a system 
with the origin left above and the y-axis directed downwards. I applied turtle graphics integrating two 
turtle commands for forwards and clockwise rotating ( vorw and rechts).“ 


The next page demonstrates Christoph Kuderli’s program and procedure. 
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Comments on (recursive) graphics programming with TI-Nspire 

Christoph Kiiderli, Switzerland 

You need the language Lua. There is an editor for LUA which is running on Nspire: Oclua (= on cal- 
culator LUA). Run Nspire and open oclua.tns and insert a new Notes page. Here you can write the 
(LUA-) program. Copy and paste the full text into the first page of oclua.tns (Ctrl C and Ctrl V). The 
output of the program - the curve - will appear immediately! 




function on.paint(gc) 
function vorw(d) 

local dx=d*math.sin(b) 
local dy=d*math.cos(b) 
local xe=xa+dx 
local ye=ya-dy 
gc:drawLine(xa,ya,xe,ye) 
xa=xe ya=ye 
end 

function rechts(winkel) 
phi=phi+winkel 

b=math.pi*phi/180 --Bogenmass 
end 

function zeichne(l) 
if l<4 then 
vorw(l) 
else 
1=1/3 

(continued in the right column) 

The oclua-Editor can be downloaded from 


zeichne(l) rechts(-60) 
zeichne(l) rechts(120) 
zeichne(l) rechts(-60) zeichne(l) 
end 
end 

xa=10 ya=150 
phi=0 b=0 

rechts(90) zeichne(300) 
end 



Tl-website: 


http://www.ticalc.org/archives/files/fileinfo/440/44075.html 
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Using TI-Nspire CAS Technology in Teaching 
Engineering Mathematics: ODEs 

Gilles Picard, Chantal Trottier, Michel Beaudin, Genevieve Savard 

Service des enseignements generaux 
Ecole de technologie superieure 
Canada 

michel.beaudin@etsmtl.ca ., gilles.picard@etsmtl.ca ., genevieve.savard@etsmtl.ca ., chantal . trottier@etsmtl . ca 

ABSTRACT 

In September 2011, our engineering school has adopted the TI-Nspire CAS CX calculator (and soft- 
ware) for all new students entering bachelor level. With both platforms and integrated facilities, 

Nspire CAS is becoming a more complete solution than Voyage 200. Although we already started to 
use this new device in the classroom, we will have a transition period where both the TI Voyage 200 
and the new CX handheld calculators can be found on the students ' desks. This talk will demonstrate 
examples of how this technology is used in classrooms when teaching differential equations. 

Nspire CX handheld with its faster CPU and extended graphing capabilities gives us even more op- 
tions for doing math in our classes. Starting with first order equations, the use of an Nspire-CAS cal- 
culator makes it easy for students to explore slope fields, numerical solutions with Euler ’s method or 
even Picard's approximations for well behaved equations. With Nspire, we can now see graphically 
how a numerical solution fits when compared to an exact solution. Consider applications of first order 
ODE's, Newton 's law of cooling for example. With Nspire, students can hook up a temperature probe 
to their calculator and gather data which will then be compared to the theoretical solution. 

Looking at second order equations, a CAS calculator can help in computing part of the algebraic so- 
lution. In the method of undetermined coefficients, we still teach the classical approach but let stu- 
dents use their handheld to substitute the candidate in the equation and solve the resulting system of 
equations. When teaching the method of variation of parameters, we show them step by step solution 
(with the help of their calculator) instead of giving them a unique formula for the solution. 

For teachers, having this technology can be very useful for grading exams. One can construct a small 
program that will yield step by step solution for the variation of parameters solution. If a student 
makes an error at the beginning (wrong solutions for the homogeneous equation), the teacher could 
easily obtain what should be the remaining solution for this student. We will demonstrate how simple 
it is to create such a program. 

Our students must learn skills in regards to technology: this is now part of the curriculum and will be 
presented. In exams or in take-home work, some problems will require the use of Nspire CAS (or Voy- 
age 200) technology. Students should also learn to validate manual results using these tools and un- 
derstand unexpected results obtained by CAS systems. We will present examples where the benefit of 
using this approach is obvious. Moreover, the majority of students will spend a lot of energy exploring 
these math concepts if these tools are used. Everything is more visual and they love to play with tech- 
nology, as long as they are told which direction to follow. 
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1. First order equations 

Consider the differential equation — = y • (4 - y) , with the initial condition y(0) = 0.2. 

dx 

We want to find the value of y when x = 1 . 

We will graph the slope field of this equation, the numerical solution (Euler method) and the exact 
solution, all in the same window. Also we will approximate the solution with a polynomial, using 
Picard’s iterations; we will graph Picard’s 3 rd iteration. 

We will also show the required value y{ 0) estimated by Euler method (in a table of values) and com- 
pare it with the exact value. 

The left screen [1] gives the slope field, and the right screen gives the slope field with the numerical 
solution, by Euler, using a plot step of 0. 1 : 




Again, the slope field and the numerical solution, but with a plot step of 0.01 : 




^Unsaved ^ 




y r ' 
\ \ \ 


\ \ \ 
\ \ \ 


\ v 
\ \ 


\ r\ \ n r 
\ \ \ \ \ \ \ 


r\\\\\\\\\\\\\\\ 

wwwwwww 



/ / / / / 

/ / / / / 

/ / / / / 

/ / / / / 

/ / / / 


/ / / / / S'/ /////// 

'////f/////J/// 
//////' '///*'// / 

/ / /V 1 - 1 /////// 
/////////// 


//////////// 


Let’s now compute the exact solution; this is possible because the differential equation is a simple 
(separable) equation and the calculator has a command for finding exact solution of well-behaved 
1 st and 2 nd order equations. 

We graph this exact solution, comparing it with the Euler solution. We come back with a plot step of 
0.1 for the numerical solution since it will help us see clearly the difference between the 2 solutions: 




Picard’s iterations 


Let’s recall that given the differential equation — = / (x,y) subject to the initial condition 

dx 


y(x o) =yo , the formula for Picard’s iterations is y n+l (x) = y 0 + J / ( t,y n (t))dt . 

*o 

In our example we have /(x, y) = y • (4 - y ), x 0 = 0 and y 0 = 0.2. 


On the left screen, we first store the function / (x, y), and then we prepare to get the iterations. On the 
right screen we press ENTER several times to iterate this formula. 




Our goal is to graph an approximate polynomial solution (see the right screen below). We chose 
Picard’s 3 rd iteration which is a 7 th degree polynomial. 




Looking at this last graph, we see that, as expected, the polynomial solution does a good job near 
x = 0 but is much in error when the x-value grows greater that 1 . 
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Table of values 

We now want to compare the numerical values obtained by these different methods. That’s why we 
use a table of values: 



The first column gives the exact (6 digits) values of the solution, so thaty(1 ) = 2.96737 ; the second 
column is the Picard’s 3 rd iteration : 2.85069, and the third column gives the value obtained by Euler, 
using a step size of 0.1 : 2.71611. 


There is another numerical method we can use, which is far better than Euler, and this is the special 
version of the Runge-Kutta method implemented in TI-Nspire (a 3 rd order adaptive step version). As 
we can see below, if we change the solution method to “Runge-Kutta” in the “Differential Equation 
Settings” window, the values in the 3 rd column of our table will change. We see below that the esti- 
mated value for y(1 ), by Runge-Kutta, is the same as the exact solution when values are rounded up to 
two decimals. When rounding up to three decimals, the difference between the two values is only 
0.001. Thus we observe a good improvement in the numerical solution and we notice that the graph 
obtained with this method is very close to the one of the exact solution! 
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2. Newton’s law of cooling, using a temperature probe 

Here is a classic example of an application of 1 st order differential equations in real life; we want to 
analyze the temperature of a liquid that is cooling down. 


Newton’s law states that this situation will be modelled by the equation — = k-(T -T E ), where T is 

dt 

the temperature of the liquid, T E is the ambient temperature (supposed constant) and k is a parameter 
depending on the environment of the experiment. The solution is T(t) = T E +C • e kt where C is the 


difference between initial temperature and ambient temperature. 


With the new TI-Nspire CX CAS, we can get real experimental data in the calculator using a tempera- 
ture-probe. Knowing the ambient temperature, we will use two of the observations to get the theoreti- 
cal solution and see how it compares to the data collected. 


Here is the experiment: We took temperatures of water cooling, from 100°C. You can see on the pic- 
ture the test tube filled with hot water, the probe connected to a calculator and a fan to keep ambient 
temperature constant. We set the calculator to take a temperature reading every 30 seconds, for 15 
minutes. 



The results are shown below. As we can see, we started a little late; water had already cooled below 
95°C. The ambient temperature was 24.0 °C. 
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On the left we have the temperature every 30 seconds. On the right, the points are some of these val- 
ues (about 1 over 4), and the software automatically connects them. 

To find these values, we solve a system of 2 equations. Let’s take data at t = 0 and 900 seconds. This 
is an arbitrary choice, we could choose other points which would give us similar, but not necessarily 
the exact same values for C and k. 

We work with time in minutes instead of seconds, so that k is not too small. We have, taking tr for 
temperature: tr{ 0) = 92.4 and tr( 15) = 35.8. 



Thus we find that tr{t ) = 24 + 68A-e °' U7l52t . This is the solution according to the theoretical model. 
We use this function to check whether it gives the results we had obtained by experiment at 300, 600 
and 30 seconds. 
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Let’s recall the temperatures obtained by 
experiment: 

71[300) = 61.9 

r(600) = 45.2 

and 

7X30) = 88.1 


We see that the results are quite satisfying, the theoretical solution giving values near what was ob- 
served. 

It is important, while we do the experiment, not to move or touch the recipient and we need to have a 
fan to insure a constant temperature for the surrounding environment. Look at the curve below. It 
represents the results of another experiment we had done earlier where we had moved the probe, and 
touched the sides of the recipient with it while having no fan. 
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3. Second order equations 

A - Undetermined coefficients 

Let’s consider the equation — ^ + 2— + 5v = 61* sin(9^). 

dt 2 dt 


We want our students to solve this equation by the method of undetermined coefficients. We show 
them how their CAS calculators can help in doing some parts of the solution. 

First they have find the roots of the characteristic polynomial and write the homogeneous solution. 



They then write the candidate for the particular solution, this candidate is A sin (9t) + B cos(9 1). Now, 
the problem is to find the values of A and B so that this candidate gives a solution when substituted in 
the differential equation. 

They can substitute it in the equation using the operator op( ) we have shown them to automate this 
substitution, and then solve the resulting system of equations. 
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Thus they obtain the general solution y g = — cos(9 1) - — sin(9 1) + • (c, • cos(2^) + c 2 • sin(2^)), and 

they even can verify that the solution is good by substituting it in the differential equation, using op( ). 


As we saw earlier, we could use a command of the calculator to get the solution in one step. The next 
screen shows this but, notice the huge difference between our previous solution and the general solu- 
tion obtained with the DeSolve command. 


Using the op( ) operator, we can verify that this long solution is indeed a solution of the differential 
equation, with a little use of trigonometric identities : 



(Comment: The DERIVE users will remember: ►sin (end of the sixth row) works similar to Trig Pow- 
ers: Sines in the Mode settings, Josef). 


B - Variation of parameters (Lagrange method) 


We ask the students to solve the differential equation 


^~Y + 6— + 9y - 45t 4 -e' 3i , 
dt dt 


using the method 


of variation of parameters. 

First they have to find the roots of the characteristic polynomial and write the homogeneous solution: 
y h = Cj t/j + c 2 u 2 . According to this method, a particular solution will be of form y p =L l u l + L 2 u 2 , 

where L x and L 2 can be variable (functional) coefficients. They have to findTi andT 2 - 


L x u l +L 2 u 2 =0 
I L[u[ + L 2 u 2 =45 t 4 -e~ 3t 


The system of equations to solve is 
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They obtain L' and , integrate to find L x and L 2 , and give the desired solution. 

We could ask them to do all this by manual calculation (classic textbook approach) but let’s get a little 
help from the CAS calculator. 



We could give them a direct formula to get this result but we prefer a step by step approach. 

We made a program with NspireCAS to test these results, even if students make a mistake, for exam- 
ple if they take the wrong homogeneous solution. This program outputs the intermediate steps we 
want to see in the students’ solutions. 
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If they make a mistake in finding the homogeneous solution, for example if they forget the middle 

term 6—, they will have y h = c x cos(3 1) + c 2 sin(3t) and their particular solution will not be good 
dt 

either, but we still can grade their solutions. 
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4. Mass-spring system 


The harmonic motion is often modelled with a mass-spring system, as we see in this picture: 



* V+ 



An object of mass m is hanging at the end 
of a spring. 

We set an x-axis such that x = 0 at rest, and 
x is positive when the spring is stretched. 


Let’s take a mass-spring system that leads to the equation 

x" + 2x + 26x = 2S(t - 3) + 2 8{t - 5), x(0) = 0, x'(0) = 1 
As we saw, x is the position of the mass so that x’ is its velocity and x” is its acceleration. And the 
Dirac Delta function S(t- a) could be visualized as hammer blow given on the object at t = a , so that 
we have two hammer blows, the first-one after three seconds, and the second-one after five seconds. 

We ask our students to solve this differential equation, with Laplace transforms. They are allowed to 
use a program they installed on their calculators. This program from Lars Fredericksen was originally 
done for the TI-89 and it was named « Laplace ». It has been modified for TI-Nspire by Philippe 
Fortin (Lycee Louis Barthou - Pau) [2] , and re-named « specfunc ». Finally the comments have been 
translated in French by Chantal Trottier, and the program received the name of « ETS-specfunc ». 

Once they get the solution, we ask them to graph it. This presents some difficulties. 

They have to get the solution, then define the function u(t ), which is the unit-step function (Heaviside 
function); its value is 0 if t < 0, and 1 if t > 0. They don’t want this function to be understood by the 
calculator when they use this program, so they will have to delete it after use. Or they could create a 
permanent Heaviside function in step(t) and replace in the output solution all occurrences of the func- 
tion u(t ) by step{t). 


ets_specfunc \soIved 

1 2 

— (x(r))+2 — -(x(/))+26 x(')-2 $(r-3)+2 sM,{xM.0, 1 } 

j,2 dt 

\dt 

"please s€ 

stmode: complex format=rectangular" 

ets_specfunc \solved 

—U‘)h2 f-M'))+26 x(f)= 2 5(/-3)+2 s(/-5).{x(/),0. 1 } 
\dt 2 dl 

/ \ 2 e J sin(5 /-15) w(/-3) 2 ( 

p • sin(5 /-25) w(/-5) sin (5 t ) w(/) 


5 e / S e f S-e / 

t | 2 e J sin(5 /— 15) w(/-3) 2 sin(5 /— 25)- i/(f-5) sin(5 /) w (*) | 

;in(5 ^-25) w(/-5) e 1 sin(5 t) u(t) 


5 ^ 5 e r 5-e / / 

2 e 3_/ sin (5 /-15) w(/-3) 2 e 5_/ s 

5 5 5 

U (/):4°'' <0 D ° ne 

lu>o 

1 
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It remains only to graph it, and to answer some questions to demonstrate their comprehension, such as 
« Is the object going down or up when the first hammer blow occurs? », « And the second? », « What 
is the greatest distance from equilibrium, when does the object reach it, and is the spring stretched or 
compressed at that point? » 



These were some examples of how we can benefit from the use of CAS technology when teaching a 
basic course in differential equations. We still teach the classic topics using a paper and pencil ap- 
proach for simple problem. We then show how technology can be used to aid in long and tedious 
computations. This gives us time to explore even more mathematics with students. 


[1] All solutions and screens were done in the original paper using TI-NspireCA CX, with O.S. 3.1 . 

I updated all screens using the latest Nspire version which is O.S. 3.4, Josef 

[2] Philippe is a very good contributor to this website in France: 
http://www.universs-ti-nspire.com/ 

You can download two respective pdf-files how to use the library (in french): 

http://seg-apps.etsmtl.ca/nspire/documents/transf%20Laplace%20prog.pdf 

https://cours.etsmtl.ca/seg/mbeaudin/MAT265/BlocLabosTI/Mat%20235%20et%20TI-Ex.supp.- 

Fonctions.pdf 

The libraries specfunc and ets_specfunc as well are among the zipped MTH96 files. 
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Spirals of Polygons 

Gunter DreeBen-Meyer & Josef Bohm 


Last September I received a mail from Ger- 
many asking how to produce spirals of poly- 
gons. Gunter DreeBen-Meyer from Berlin com- 
plained that this was very easy done with DE- 
RIVE applying one single VECTOR-command. 

Students could apply matrix calculations first 
and then individually proceed with the so gen- 
erated graphs (colours, ...). But unfortunately 
he and hiscolleagues were unable to create the 
spirals with TI-NspireCAS. 

His mail ended with: 

Do you have any idea or advice how to achieve 
the spirals with Tl-Nspire? 



Best regards, A Spiral of (DERIVE-) Squares 

Gunter DreeBen-Meyer, Berlin 


I wrote back that I would try to find a solution. I must admit that I didn’t remember at this time 
that Maria Koth had treated these spirals in DNL#33. So I started from the very beginning 
and my first solution might be interesting but it is not very elegant at all. The problem I en- 
countered was - again - missing the plotting commands of Voyage 200. The only one way 
plotting segments is via one - or even more - scatter diagrams ... 

I sent file Poly_Spiralen.tns: parameters must be entered on a Notes page. 
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Enter the parameters on this Notes page: 
starting point pt, number of vertices of the regular polygon n and angle of rotation a 

pt:=[() 2] ► [0 2], 11 ► 5. , a ► 0.314159 

We can install sliders for n and for a. 

That's all, don't enter any other data, calculation will be done in the background! 

xw:=poly_x(pt,n) ► { 0. ,-1.9021 1,-1. 17557, 1.17557,1. 9021 1, 0. } 
yw:=poly_y(pt,n) ► { 2., 0.618034, 1.61803, 1.61803, 0.618034, 2. } 
fa:=q(n,a) ► 0.850651 


xw and yw are the vertices of the first polygon. The coordinates of the next 10 generations 
are given in the - flexible - spreadsheet page - up to a 20-gon (xwl ,yw1 ) to (xwl 0,yw1 0). 
Then 1 1 scatter diagrams are defined - one for each generation. 


I print the handheld view to save space. 


► *Poly_Spiralen ^ <HP 


< j > 1 1» =S. I < m a =14 

scatter plots/ 




<1 

HiniTi ► *Poly 

_Spiralen ^ 

<IO 

♦ 

A 

B 

c xwl 

D ywl 

|0 

1 

= 

= XW 

=yw 



13 

0.929446 

1 .77091 

0.525496 

1 .87048 

0 

14 

0. 

2. 

-0.40395 

1 .90044 

15 





16 












1 =1 fa (d7 7^ sin(a)+£?7^- cos(o)) 



This is in fact not very comfortable - compared with the DERIVE-VECTOR procedure. 
I will try a LUA-script - but this will last a little while. 
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From Berlin: 

Your proposal looks quite interesting, in particular storing the vertices in the spreadsheet. 
I will inspect it more accurately. 

Many thanks and best regards 

Gunter DreeBen-Meyer 


Dear Mr. DreeBen-Meyer 

This is the first solution of your problem worked out with the script language LUA. It was sent 
by a Swiss DUG-member, Alfred Roulier. You can see how excellent cooperation is working 
in the DUG-community, btw, many thanks Alfred. 


I am working on a LUA-script calculating and plotting the single polygons stepwise instead of 
generating the whole point matrix - which was Alfred’s way. 


This is Alfred Roulier’s file whirl. tns (handheld screen shots). The important part of the LUA- 
script is given on the next page. Program whirl() generates the numeric values which are 
then used by LUA. 


1.1 | 1 2 1 1 3 

► whirls 


1 i-i 1 

1 12 | 1.3 


ne:=4 ► 4 Anzahl Ecken -number of vertices 
ns: = 12 ► 12 Anzahl Schritte -number of steps 
si: = 100 ► 100 Seitenlange -side length 

dip: =0.1 ► 0.1 
whirlQ ► Done 
r 

► { 100,91.3377,83.4258,76.1993,69.5987,6 


0 




C\ A 1 A 1 A 


1 .1 A L A r-' A 7 A S A Cl 1 / M 


whirl ^ 


"whirl" stored successfully 


Define whirl()= 

Prgm 
Local ij 

xk. =newMat {ne f ns):yk. =newMat(tf 
*: =ne w Lis t{ns) : ?{ 1 ] : =sl 
^:=seq(r d<pj,0 t ns-l) 

© Radien berechnen- Calculation of radii 
For i, l,ns-l 

n 




cos 


0 


!► whirls <HQ 


5 

1.2 1.3 1 

Fl 



The vertices are lying on a logarithmic spiral with the centre of the polygon as pole, 
(from A. Roulier, whirl. doxc inspired by an article in Scientific American, 1965) 
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function on.paint(gc) 

ne=(var.recall("ne") or 1) — import number of vertices 
ns=(var.recall("ns") or 1) — import number of steps 
— Punktmatrizen definieren 
xk = {> ; yk = {> 
for j=l,ne do 

xk[j] = {> ; yk[j] = {} — create a new row 
for i=l,ns do 
xk[j][i] = 0 ; yk[j][i] = 0 
end 
end 

— import point coordinates 

xk=(var.recall("xk" or 2)) ; yk=(var.recall("yk" or 2)) 

— plot polygons 
for i=l,ns do 

for j=l,ne-l do 

gc:setColorRGB(255*(l-(j-l)/ne),0,255*(j-l)/ne) 
gc:drawLine(w/2+xk[j][i],h/2-yk[j][i],w/2+xk[j + l][i],h/2-yk[j + l][i]) 
end 

gc:setCoiorRGB(0, 0,255) 

gc:drawLine(w/2+xk[ne][i],h/2-yk[ne][i],w/2+xk[l][i],h/2-yk[l][i]) 

end 

end 


Again a mail from Berlin: 

Dear Mr Bohm, 

many thanks for the mail and the files. 

How easy was it with DERIVE! And the issue mappings and matrices was reproducible and 
investigable by students. It’s a pity that these fine properties have not been taken over to 
NSpire! 

I should start programming. Many things seem to become possible. I see that you have been 
grasped by programming, too. 

Regards from the North Sea, 

Yours G. D.-M. 


Dear Mr. DreelJen-Meyer, 

finally it worked. I was able to transfer my iterative algorithm into a LUA-program (see the 
attached file Poly_Spirals_LUA.tns) 

Two nice problems for the students could be offered. 

How to generalize the mapping that 

(1) the direction of the spirals are reversed, and 

(2) the spirals are not inscribed, but circumscribed? (possible answers see page 27) 
Best regards, Josef B. 


p 24 


G. DreeBen-Meyer & J. Bohm: Spirals of Polygons 


DNL 96 


Die Eingabe der Parameter n, a und its erfolgt uber die Schieberegler auf page 3 
Enter the parameters n, a and its on page 3 using the sliders 

Startpunkt - initial point 

pt=[o 3] • [0 3] 

n ' 8. 
its ► 6. 

ct - 47. 

Nichts mehr eingeben - Nothing to enter! 

xw:=raund(l000- (so- poly_x(ptn))^)/1000 ► { 0. ,-106. 066,-150. ,-106. 066.0. ,106. 066, 150.406.066,0. } 
yw: =round(lQ00 ■ (so ■ poly_y(ptn)),3)/1000 ► { 150., 106. 066,0. ,-106. 066,-150. ,-106. 066,0., 106.066,150. } 

Daten fur das Script (nichts eingeben) - Data for the script, nothing to enter: 

a:=approK[ j ► 0.820305, fa:=apprax(q(n,a)) ► 1.0153 

\l80/ 
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This is the complete LUA-script: 

platform. apilevel='2. 2' 
function on.construction() 

h=platform.window:height(); w= platform. window: width() 
end 

function on. resize(width, height) 
h = height; w=width 
end 

function on.enterKey() 
platform, window: in validate() 
end 

function on.paint(gc) 
n=var.recall("n");a=var.recall("a") 
xw = {}; yw={> 
for i=l,n+l do 
xw[i]=0;yw[i]=0 
end 

xw = va r. reca 1 1 ("xw") ; y w= va r. reca 1 1 ( "y w" ) 
for i=l,n do 

gc: setColorRGB(0, 0,255) 

gc: drawLine(w/2+xw[i],h/2-yw[i],w/2+xw[i+l],h/2-yw[i+l]) 
end 

its=var.recall("its");fa=var.recall("fa") 

xn={};yn={} 

for i=l,n do 

xn[i]=0;yn[i]=0 

end 

— the iterative algorithm 
for j=l,its do 

for i=l,n+l do 

xn[i]=fa*(xw[i]*math.cos(a)-yw[i]*math.sin(a)) 

yn[i]=fa*(xw[i]*math.sin(a)+yw[i]*math.cos(a)) 

end 

for i=l,n do 

gc: setColorRGB(0, 0,255) 

gc: drawLine(w/2+xn[i],h/2-yn[i],w/2+xn[i+l],h/2-yn[i+l]) 
xw[i]=xn[i];yw[i]=yn[i] 
end 

xw[n+l]=xn[n+l];yw[n+l]=yn[n-i-l] 

end 

— end of iteration 
end 

Exchanging the gc:setColorRGB(0, 0,255) by gc: setColorRGB(10*j,100+5*j, 255-1 0*j) results 
in a nice shading of the nested polygons. 
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The final solution - without LUA came in from Berlin: 

... many thanks for your files and your intellectual efforts. I will study the program code, 
and some days later 

Hello Mr. Bohm and Karl-Heinz Keunecke, 

I have solved our problem (matrices and spirals) in another way. 

The vertices of a polygon are defined by polar coordinates. The coordinates of the initial po- 
lygon are stored together with the coordinates of the rotated figures in one list for the 
x-values and a second one for the y-values. 

The scatter plot of the lists gives the plot. The attributes must be adopted. Calculation of both 
lists is performed within the Nspire environment (not in a LUA-program). 

Best regards from Berlin 
Gunter D. 


I intend to write a short article for the Tl-News. 


► *Document2 ^ 


"programmzeilen" stored successfully 


Define programmzeilen()= 

Prgm 

Local ij.xw:={ i_i i.i } 

For 1, l,dz 

For i, l,/i+l | 

^ x>p[z+(/-l)- (/?+l)]:=cos(/wvr 7+(/-l)- dw)- r 

yu[i+{l- 1)- (tf+l)]:=sin(mv- ?-»-(/— l) dw)- r ] 
EndFor: EndFor 
End Prgm 


E3 


0 


] 11 | 1 .2 | 1.3 | 


► *Document2 


Anzahl Ecken n:=5 * 5 Diehungen dz:=12 * 12 

K 

Drehwinkel d w :=1 5 * » — — Radius r :=6 *6 
12 

2 k 2-tc 

* — Mittelpunktswinkel 

n 5 

Beiechnung 




>/r+i 


S (H 4C0S (ir) 


progtammzeilenQ * Done |*usfiihren des Progiamms 


0 


0 



<0* 



This is the family of pentagons with dw = 80°. 
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Gunter D.-M. derives the formula for the 
factor fin the same way as Maria Koth did 
(DNL#33). He shows a pretty picture 
accomplished by one of his students adding 
a nice pattern to a spiral of squares. 



As this is still the DERIVE newsletter I’d like to add some DERIVE made plots, too. On page 
23 I posed possible problems for the students. 

d 


P0L_R0T(r , n, Q, j, d 1) ITERATES (S (n , 



inscribed 


POL_ROT 3 , 5, , 5, -1 

l 18 



<t0 ■ ROT (cp , p), p, POLYCr, n), j) 




circumscribed 


circumscribed 
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In the last DNL I recommended Benno’s wonderful book - in German only and until 
yesterday for MS-Excel. I incited Benno to rewrite this great collection of computer 
simulations for Tl-Nspire - and just today a mail came in: 


“I finished your challenge just now. There are 92 Nspire files - many of the including LUA- 
scripts. I am sending the preliminary text as a pdf-file.” 
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This is the flyer for the Excel-edition. The next page shows some graphs 
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I present some screen shots from Chapter 18 “Tennis”. There are 24 chapters on nearly 200 
pages. 





siege 



- 

nrm» =.55 1 <n>n i =197. 

Wins of A in %: 

1 

Number of Games : 2 5 

64 

1 


v 

0 


You can see the probability tree for the run of 
one game. 

The last row shows the outcome, either on 
the left (SA = A wins the game) or on the 
right (SB = B wins the game). 

The positions in the last row are: 

Win A, Advantage A, Deuce, Advantage B 
and Win B. 


The respective Excel-Graph offers more space for annotating the knots in the graph. 
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Reste und ihre Muster 
Complex Numbers on Norwegian Sweaters 

Heinrich Ludwig, Germany 


Die graphische Darstellung mathematischer Prozesse und Strukturen fuhrt oft zu asthetisch 
ansprechenden Bildern. Fraktale Mengen wie Julia- Oder Mandelbrotmengen beispielsweise 
fanden ihren Weg in die popularwissenschaftliche Literatur hauptsachlich wegen ihrer reiz- 
vollen Abbildungen. Aber auch die elementare Zahlentheorie kann interessante Grafiken 
erzeugen. Derfolgende Text mochte das mit einem Beispiel belegen. Einen kleinen Vorge- 
schmack geben die abgebildeten Grafiken. 


Graphic representation of mathematical processes and structures often leads to aesthetic 
appealing pictures. Fractal sets like Julia- or Mandelbrot sets found their way into popular 
literature primarily because of its attractive pictures. But even elementary number theory is 
able to produce interesting graphics. The following article shall demonstrate this. Take the 
pictures below as a foretaste. 



'V* 
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Reste in Z und Z[i], Remainders in Z and Z[i] 

Die Menge der ganzen Zahlen Z ist ein euklidischer Ring. Damit ist u.a. gesagt, dass es in Z 
eine Division mit Rest gibt. Legt man einen bestimmten Divisor (= Modul) m fest, dann ist 
f(x) = MODS(x,m) eine Funktion, die jeder ganzen Zahl ihren Divisionsrest zuordnet. Dabei 
ist der Betrag des Funktionswerts, d.h. des Restes, kleiner als der Betrag des Moduls m. 
Diese Eigenschaft ist fur das Weitere wichtig, denn sie stellt sicher, dass f(x) nur endlich vie- 
le Funktionswerte hat, obwohl f fur unendlich viele x-Werte definiert ist. 

Auch die Menge aller ganzen komplexen Zahlen Z[i] ist ein euklidischer Ring, d.h. die von 
den ganzen Zahlen Z her bekannte Division mit Rest gibt es auch in Z[i]. Leider arbeitet die 
von DERIVE bereitgestellte MODS-Funktion nicht mit komplexen Zahlen. Das macht nichts; 
eine Funktion, die den komplexen Rest berechnet, kann man mit DERIVE programmieren. 

The set of the integer numbers is an Euclidean ring. Among other properties exists a division 
with remainder in Z. Fixing a certain divisor (= modulus) m, then f(x) = MODS(x,m) is a func- 
tion assigning every integer x its remainder after division by m. The absolute value of the 
remainder (function value of MODS) is smaller than the absolute value of the modulus. This 
is important for the following because it makes sure that f(x) has only a finite number of func- 
tion values although f is defined for infinite many x-values. 
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The set of all integer complex numbers Z[i] form an Euclidean ring, too, i.e. the division with 
remainder is also existing. Unfortunately function MODS provided by DERIVE does not sup- 
port complex numbers. Doesn’t matter, a function which calculates the complex remainder 
can be defined by DERIVE. 


M0DC(a, m, c, d) := 
Prog 

If m = 0 


RETURN a 
c := a/m 

d := FL00R(RE(c)) + L ■ FLOOR (IM(c)) 
d := a - m ■ d 
If RECd) = 0 
d := IM(d)-L 
If IM(d) = 0 
d := RECd) 

RETURN d 


MO DCS (a, m, c, d) := 
Prog 

If m = 0 


RETURN a 
c := a/m 

d := FLOOR (RE (c) + 1/2) + i . FL00R(IM(c) + 1/2) 
d := a - m-d 
If RECd) = 0 
d := IM(d)-L 
If IM(d) = 0 
d := RE(d) 

RETURN d 


Die beiden IF-Anweisungen am Ende von MODC bzw. MODCS sind aus mathematischer Sicht uber- 
flussig. Sie verhindern nur die Auswirkung eines internen Fehlers von DERIVE 6. 

Both IF-commands at the end of MODC and MODCS are not necessary from the mathematical point 
of view. They prevent an internal error of DERIVE. 

MODCS(a,m) und MODC(a,m) sind die komple- |viqdc(5 + 3.l 4 + L) - l + 2 -i. 

xen Analoga zu MODS und MOD, wobei a der 

Dividend und m der Modul ist. Ein Beispiel: M0DCS(5 + 3 -l, 4+l)=1+2*l 


MODCS(a,m) and MODC(a,m) are the complex 
analogues of MODS and MOD, where a is the 
dividend and m the modulus. An example: 


M0D(5 + 3-l, 4 + l) = 1 + 2-1 
MODS (5 + 3-i, 4 + l) = 2 — 2-L 


Comment: See the definitions of MOD and MODS given in the DERIVE Online-Help: 

MOD(m, n) simplifies to m modulo n (i.e. the nonnegative remainder of m/n) . 

MODS(m, n) simplifies to the symmetric m modulo n which will be in the half-open inter- 
val [-n/2,n/2). 

Ein Modul m aus Z \ {0} sei festlegt. Wenn man die Funktion MODS mit alien moglichen gan- 
zen Zahlen futtert, dann bekommt man als Wertemenge alle ganzen Zahlen des Intervalls 
(- m/2, m/2], Beispielsweise bekommt man fur m = 13 folgende Restemenge: {-6, -5, -4, -3, 
-2, -1,0, 1,2, 3,4, 5, 6} 

A modulus from Z \{0} be fixed. Feeding function MODS with all possible integers then one 
receives all integers e (- m/2, m/2], E.g. for m = 13 the set of remainders is given by {-6, -5, 
-4, -3, -2, -1,0, 1,2, 3,4, 5, 6}. 

VECT0R(M0DS(k , 13), k, 100, 120) 

[-4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3] 

VECT0R(M0D(k, 13), k, 100, 120) 

[9, 10, 11, 12, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 0, 1, 2, 3] 

Was erhalt man bei analogem Vorgehen im Ring Z[i]? Ein Modul m aus Z[i] \ {0} sei festge- 
legt. Futtert man die Funktion MODCS(x, m) mit alien moglichen komplexen ganzen Zahlen 
x, so treten auch hier nur endlich viele Reste auf. Der Grund dafur ist derselbe wie in Z. 

Fur m = 4+i und eine Auswahl von Dividenden x = a+b-i mit -5 < a < 5, -5 < b < 5 sind die 
Ergebnisse in folgenden Tabelle (Matrix) zusammengefasst. 
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What is happening in ring Z[i]? A modulus m e Z[i] \ {0} is fixed. Feeding function MODCS 
with all possible complex integer numbers we receive again a finite number of remainders. 
The reason is the same like in Z. See m = 4 + i and a selection of dividends x = a + b i with -5 
< a < 5, -5 < b < 5. Then the result can be given in form of a table of remainders: 


VECTOR (VECTOR (MODCS (a + b-L, 

4 + l)j 3l, 

-5, 5), 

b , 5, -5, 

- 1 ) 




2-L 1 + 2-L 

-2 + L 

-1 + L 

L 

1 + L 

-2 -1 

0 

1 

2 

L 1 + L 

-2 

-1 

0 

1 

2 -1 - L 

— L 

1 - L 

2 - L 

0 1 

2 

-1 - L 

—L 

1 - L 

2-L -1 - 2-L 

- 2-L 

2 -l 

1 + 2-L 

— L 1 - L 

2 - L 

-1 - 2-l 

- 2-L 

2-L 

1+2 -L -2 + L 

-1 + L 

L 

1 + L 

- 2-L 2 -l 

1 + 2 -l 

-2 + l 

-1 + L 

L 

1 + L -2 

-1 

0 

1 

-1 + L L 

1 + L 

-2 

-1 

0 

1 2 

-1 - L 

—L 

1 - L 

-1 0 

1 

2 

-1 - L 

—L 

1 - L 2-L 

-1 - 2-L 

- 2-L 

2-L 

-1 - L — L 

1 - L 

2 - L 

-1 - 2-l 

- 2-L 

2-L 1 + 2-l 

-2 + L 

-1 + L 

L 

-1-2 -L - 2-l 

2-L 

1 + 2-L 

-2 + L 

-1 + L 

l 1 + L 

-2 

-1 

0 

-2 + L -1 + L 

L 

1 + L 

-2 

-1 

0 1 

2 

-1 - L 

—L 

-2 -1 

0 

1 

2 

-1 - L 

-L 1 - L 

2 - L 

-1 - 2-L 

- 2-L 

Comment: Let’s compare the first rows of this matrix created using DERIVE’s MODS-function: 

VECTOR(VECTOR(MODS(a + b - L , 4 + L ), a , - 5 , 5 ), b , 5 , - 5 , - 1 ) 

' 1-2 -L 2 - 2-l 

-1 - 3 - i . 

-3-L 

1 - 3* L 

2 - 3-L 

— 1 — 4-L -4-L 

0 

1 

2 

1 - 3-l 2 - 3 ■ l 

-1 - 4-l 

- 4-L 

0 

1 

2 -1 - L 

—L 

1 - L 

2 - L 

0 1 

2 

-1 - L 

—L 

1 - L 

2-L -1 - 2-L 

- 2-L 

1 - 2-L 

2 - 2-l 

— L 1 - L 

2 - L 

-1-2 -L 

- 2-L 

1 - 2-l 

2 - 2-L -1 - 3-l 

- 3-L 

1 - 3-l 

2 - 3-l 

- 2-l 1 - 2-l 

2 - 2-l 

-1 - 3-L 

- 3-l 

1 - 3-l 

2 - 3-L -1-4 -L 

MODS (-5 + 5-L, 

- 4-l 
4 + L) = 

0 

1 - 2-L 

1 


Comment: So we have different symmetric re- 
mainders. I checked by subtracting both remain- 
ders from the dividend and then calculating the 
MOD-, MODS-, MODC-, and MODCS-values. 

Comment: I am sorry but I cannot give an expla- 
nation for this. I am sure that we have several 
complex number specialists among us who can. 


MO DCS (-5 + 5-l, 4 + l) = Z*l 
M0D(-5 + 5-l - 1 + 2-l, 4 + i.)=0 
M0D(-5 +3-i., 4 +l) = 0 
MO DCS (-5 + 5-l - 1 + 2-l, 4 + 0= 0 
MODS (-5 + 5-l - 1 + 2-l, 4 + 0= 0 
MODS (-5 + 5-l - 2-l, 4 + 0=0 
MO DCS (-5 + 5-l - 2-l, 4 + 0= 0 


Eines lasst sich sofort erkennen: Die Funktion f(x) = MODCS(x, m) ist „2-dim-periodisch“: 
Geht man in der Matrix 4 Spalten nach rechts und eine Zeile nach oben, dann erhalt man 
ungeachtet des Ausgangspunkts dieselbe Zahl. Das gleiche sieht man, wenn man eine Spal- 
te nach links und 4 Zeilen nach oben geht. Eine weitere Beobachtung macht man, wenn man 
alle vorkommenden Funktionswerte zusammentragt: Sie liegen in der komplexen Ebene in 
einem Quadrat, dessen eine Seite m ist und dessen andere m i. Fur m = 4+i findet man 
daher insgesamt |m| 2 = 17 verschiedene Werte. (Diese Feststellungen sind unmittelbare 
Auswirkungen der Funktionsvorschrift von f(x).) 
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What we can recognize immediately is that function f(x) = MODCS(x,m) is “2-dim-periodical”: 
Starting at any element of the matrix and stepping 4 columns to the right and then one row 
up we will end at the same number. The same happens making 4 steps to the left followed 
by one step down. We can make a second observation when collecting all appearing function 
values. They are within a square in the complex plane with one side m and the other one m i. 
Therefore for m = 4 + i there are totally |m| 2 = 17 different values. (These facts are conse- 
quences of function rule of f(x).) 



+ 




MODCS generated remainders and MODS generated remainders 


4 


4 ■ ■ ■ 


2 


2 ■ ■ ■ 


2 


4 


2 


4 


MOD and MODC give the same 17 points 

Quadratische Reste, Quadratic Remainders 

Interessant wird es erst, wenn man f(x) = MODCS(x, m) mit den Quadraten aller ganzen 
Zahlen futtert. Jetzt treten nur bestimmte Funktionswerte (= Reste) auf. Wahlt man zum Bei- 
spiel m = 13, so treten von den 13 moglichen Resten nur 7 verschiedene auf, namlich: 

{-4, -3, -1,0, 1,3, 4}. 

(Die Anzahl dieser so genannten quadratischen Reste in Abhangigkeit vom Modul ist eine 
sehr interessante Funktion; sie soil aber hier nicht weiter untersucht werden.) 

Die Idee, nur die Reste von Quadratzahlen zu betrachten, soil auf Z[i] ubertragen werden. 
Wir erzeugen die Matrix von oben nochmals, diesmal aber mit dem Argument (a + b i) 2 : 


It becomes really interesting if we feed f(x) = MODCS(x, m) with the squares of all integers. 
Now only certain function values (= remainders) appear. Choosing e.g. m = 13 instead of 13 
possible remainders we find only seven ones, which are: 

{-4, -3, -1,0, 1,3, 4}. 

(The number of these so called quadratic remainders depending on the modulus is a very 
interesting function, which cannot be investigated here.) 
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The idea to only observe the remainders of squares shall be transferred on Z[i], We evaluate 
once more the matrix from above exchanging a + bi by (a + bi) 2 : 

2 


(TOR (VECTOR (MODCSCCa 

+ b ■ l) 

, 4 + 

Op a, 

l n 

l n 
1 

b, 5 , 

- 5 , -1) 



L 

-2 

2 

- 2-i 

-1 

2-i 

— i 

1 

0 

1 

—i 

-1 

2-i 

-i 

1 

0 

1 

-i 

2-i 

-1 

- 2-i 

2 

0 

1 

-i 

2-i 

-1 

- 2-i 

2 

-2 

i 

i 

-2 

-1 

- 2-i 

2 

-2 

i 

i 

-2 

2 

- 2-i 

-1 

2-i 

L 

i 

-2 

2 

- 2-i 

-1 

2-i 

—i 

1 

0 

1 

- 2-i 

-1 

2-i 

-i 

1 

0 

1 

—i 

2-i 

-1 

- 2-i 

1 

0 

1 

-i 

2-i 

-1 

- 2-i 

2 

-2 

i 

i 

2-i 

-1 

- 2-i 

2 

-2 

i 

i 

-2 

2 

- 2-i 

-1 

-2 

i 

i 

-2 

2 

- 2-i 

-1 

2-i 

— L 

1 

0 

2 

- 2-i 

-1 

2-i 

-i 

1 

0 

1 

— i 

2-i 

-1 

—i 

1 

0 

1 

-i 

2-i 

-1 

- 2-i 

2 

-2 

L 


An der 2-dim-Periodizitat hat sich, wie zu erwarten war, nichts geandert. Die Anzahl ver- 
schiedener Funktionswerte hat gegenuber vorher aber abgenommen. Von den urspriinglich 
17 sind nur diese 9 geblieben: 

{-2, -1 , -2-i, -i, 0, i, 2 i, 1,2}. 

Die Wertemenge der Funktion g(x) = MODCS(x 2 , m) ist also nur eine Teilmenge der Werte- 
menge von f(x). Wahrend f(x) ein liickenloses Punktegitter innerhalb eines Quadrats erzeugt 
(s. Abb auf Seite 34), erzeugt g(x) eine Punktmenge mit etlichen „Lochern“. Ein Muster ent- 
steht (siehe die nachste Abbildung). 

2-dim-periodicity did not - as expected - change. But the number of different function values 
decreased. Instead of 17 from before only 9 of them remained, namely: 

{-2, -1 , -2-i, -i, 0, i, 2-i, 1,2}. 

The range of g(x) = MODCS(x 2 , m) is only a subset of the range of f(x) = MODCS(x,m). 
While f(x) produces a grid of points without gaps within a square g(x) results in a point set 
with several “holes", a pattern appears (see the graph). 
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Bilderzeugung, Generating the Patterns 


Zwei Dinge fehlen noch. Zum einen ist das eben erzeugte „Muster“ noch nicht wirklich inte- 
ressant. Doch das liegt nur an der Wahl des Moduls m, der in dem relativ kleinen Quadrat 
eben nur ein einfaches Kreuz entstehen lasst. Zum anderen sollte das Sammeln der ver- 
schiedenen Teilungsreste aus Matrizen wie oben gegeben einem DERIVE-Programm uber- 
tragen werden, da eine manuelle Bearbeitung muhsam ware. 


Die Aufgabe des Restesammelns iibernimmt die Funktion ShowAIIRemainders(m). Setzt 
man fur den Parameter m eine Zahl aus Z[i] \ {0} ein, so gibt die Funktion einen Vektor zu- 
riick, der (ohne Duplikate) alle Reste bezuglich m enthalt. 


OBACHT: Fur m wird vorausgesetzt, dass RE(m) > 0 und IM(m) > 0. Das bedeutet keine 
Beschrankung der Allgemeinheit, da m, m i, -m und -m i dieselben Reste erzeugen und min- 
destens eine dieser vier Zahlen die Bedingung erfullt. 


Two things are still missing. Firstly, the “pattern" presented above is not really breath taking. 
The reason for this is only the choice of the modulus m, which can create only a small cross 
within the small square - this can be changed easily. And secondly we should find a DERIVE 
program for collecting the different remainders given in the matrices above. 

The task of collecting the remainders is done by function ShowAIIRemainders(m). For any 
meZ[i] \ {0} this function gives back a vector (list) containing all remainders with respect to m 
(without duplicates). 

TAKE CARE: RE(m) > 0 and IM(m) > 0. This no restriction of generality because m, m-i, -m 
and -m-i produce the same remainders and at least one of these numbers fulfill the condition. 
2 

g(x) := x 

ShowAl 1 Remainders (m, k, j, r, s:=[0]):= 

Prog 
k := 0 
Loop 

If k = RE(m) exit 

j := 0 

Loop 

If j = RE(m) exit 
r := M0DCS(g(k + j-L), m) 

If - MEMBER?(r , s) 
s := ADJ0IN(r, s) 
j :+ 1 
k :+ 1 
k := 0 
Loop 

If k = IM(m) exit 

j := 0 

Loop 

If j = IM(m) exit 
r M0DCS(g(k + RE(m) + j-L), m) 

If - MEMBER?(r , s) 
s ADJOIN (r , s) 
j :+ 1 
k :+ 1 
RETURN s 

ShowAl 1 Remainders (4 + l) = [-l, 2, -2, 2-i., 1, -2 ■£., i, -1, 0] 

DIM(ShowAl 1 Remainders (40 + !!■£.)) - 861 
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Die Struktur der Funktion ShowAIIRemainders(m) ist leicht zu erkennen. Allerdings macht sie Ge- 
brauch von einer Umsortierung, die kommentiert werden sollte. 

Fur jeden gegebenen Modul m ist die Menge der Werte der Funktion MODCS(x,m) endlich, weil 
MODCS(x,m) fur festes m hinsichtlich x 2-dim-periodisch ist. Darum geniigt es, eine Periode zu 
durchmustern. Die Periode hat in der GauBschen Ebene die Form eines Quadrats mit den Ecken 
m/2*[-1-i, 1-i, 1+i, -1+i], (Dieses Quadrat ist in den beiden Abbildungen strichliert eingetragen.) Das 
Quadrat umfasst ||m|| Punkte aus Z[i], Weil die Seiten des Quadrats i.a. nicht achsenparallel liegen, 
eignet es sich schlecht fur eine Durchmusterung. Darum wird es nach Art des Beweises des Hohen- 
satzes in zwei neue Quadrate umgewandelt. Betrachte dazu die nachste Illustration. Die zwei neuen 
Quadrate enthalten dieselben Reste wie das ursprungliche, liegen aber achsenparallel. Den beiden 
Quadraten entsprechend besteht das Programm aus zwei Teilen, deren jeder jeweils eines der beiden 
Quadrate durchmustert. 

Structure of function ShowAIIRemainders(m) can be recognized easily. But it applies a re- 
sorting which shall be commented. 

For a given modulus m the set of function values of MODCS(x,m) is finite, because 
MODCS(x,m) for a fixed m with respect to x is 2-dim-periodical. So it is sufficient to observe 
only one period. The period has the form of a square with vertices m/2*[-1-i, 1-i, 1+i, -1+i] in 
the GaulJ plane. (This square is added in the graphs above with dashed lines.) This square 
contains ||m|| points of Z[i], Because the sides of this square usually are not parallel to the 
axes it is not very suitable for searching for points. Therefore it will be - like in the proof of 
the altitude theorem - transformed into two new squares. The graph below is illustrating the 
procedure. Both new squares contain the same remainders as the original ones but are lying 
parallel to the axes. The program consists of two parts each of them investigating one of the 
two squares. 



Die komplexen Zahlen, die ShowAIIRemainders(m) liefert, miissen noch in ein kartesisches 
Koordinatensystem eingetragen werden. Dazu muss jede Zahl in einen 2-dimensionalen 
Vektor umgewandelt werden. Das besorgt die Funktion Pattern(m). Wie ihr Name verrat, 
erzeugt Pattern(m) endlich das, was eingangs versprochen wurde: schone Muster. Ihr Reiz 
liegt darin, dass man fur (beinahe) jede Zahl m ein neues, andersartiges Muster bekommt 
und sich nur schwer vorhersagen lasst, welche Werte m besonders schone Muster ergeben 
werden. Das ladt ein zum Probieren und bringt immer wieder Uberraschungen mit sich. 

Pattern5(m) und Pattern9(m) sind zwei Erweiterungen von Pattern(m), die das Muster zu den 
4 Seiten bzw. rundherum mit Duplikaten fortsetzen, so dass es ausdrucksvoller zur Geltung 
kommt. 
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The complex numbers delivered by ShowAIIRemainders(m) must be entered in a Cartesian 
system of coordinates. Every number must be transformed into a 2-dimensional vector. This 
is performed by function Pattern(m). According to its name Pattern(m) produces pretty pat- 
terns. Its special attraction is that nearly every m results in a new and different pattern and it 
cannot be predicted in advance which values of m will give extra pretty patterns. This invites 
for experimenting. 

Pattern5(m) and Pattern9(m) are two extensions of Pattern(m), which append the pattern in 
four and eight directions to make it much more impressive. 


Pattern(m) := VECTOR ([ RE (u ) , 

f 


Patterns (m) := APPEND 


VECTOR 


Pattern9(m) := APPEND 


VECTOR 


I M C li ) ] p u, ShowAl 1 Remainders (m)) 



RE(li) 

IM(lj) 

■’i 



RE(u + m) 

IM(u + m) 




RE(u - m) 

IM(u - m) 

P li ? ShowAl 1 Remainders (m) 



RE(li + m - l) 

IM(li + m- l) 



L 

. RE(li - m- l) 

IM(li - m- l) . 

j 

j 


RE(li) 

IM(li) 



RE(li + m) 

IM(lj + m) 



RE C li - m) 

IM(li - m) 



RE(lj + m« l) 

IM(li + m - l) 



RE(u - m- l) 

IM(li - m- l) 

, Li, ShowAl 1 Remainders (m) 


RE fu + m - (l + 1)) 

IM(li + m- (l + 1)) 



RE(lj + m« (i. - 1)) 

IM(li + m- (l - 1)) 



RE(u + +1)) 

IM(li + m- £-l + 1)) 



RE fu + m-C-L - 1)) 

IM(u + m- (-L - 1)) . 

J 

j 


Wer die Ausfiihrungen bis hierher verfolgt hat, dem seien die versprochenen Muster nicht 
langer vorenthalten! Man offne ein 2D-Graphikfenster, stelle einen Zeichenbereich der Breite 
und Hohe 100 mit Zentrum bei [0,0] ein und setze das Verzerrungsverhaltnis auf 1:1. Jetzt 
sollte DERIVE noch angewiesen werden, dicke Oder mittelgrofte Punkte zu zeichnen und 
diese nicht zu verbinden. Dann lasse man den folgenden Ausdruck plotten: 

If you have followed the article until now you should be presented the promised patterns! 
Open a 2D-plot window and set the aspect ratio to 1:1. Take Point Size Large (or medium) 
and don’t connect the points. Then plot the following expressions: 

Patterns (49) 

Pattern9(41) 

Pattern9(7 + 27- i.) 

Patterns (29 + l) 

Pattern(29 + l) 
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Durch Variation des Moduls m lassen sich unzahlige weitere Muster gewinnen. Letztlich sind 
der Vielfalt der Muster nur durch die Rechenzeit Grenzen gesetzt. 

By varying the modulus m we can get numerous other patterns. Restrictions are given only 
by calculation time. 

Zum Schluss fuge ich noch einige Kommentare und Fragen an: 

Anmerkungl: 

Interessante Muster findet man insbesonders fur Moduln (=Divisoren) m, fur die gilt: 

a) m aus Z und ungerade, besonders hiibsch sind m = 7, 9, 27, 41, 49, 73. 

b) ||m|| = m m* ist Quadratzahl, (d.h. Real- und Imaginarteil gehoren zu einem pythagorai- 
schen Tripel), z.B. m = 40+9i 

Nett sind auch: m = 9+i, 15+i, 39+i, 39+5i, 39+22i, 39+29i, 40+1 1 i, 40+23i, 40+29i, 40+39i, 
41+i, 41 +9i, 41 +20i, 41+38i. 

Die Muster fur a + b i und b + a-i sind bis auf Spiegelung gleich. 

Fur a aus Z und a <> 0 (mod 4) scheint zu gelten: Die Muster bezuglich a und a+a-i sind 
gleich. 

Anmerkung2: 

Uber Hinweise auf Moduln m, die besonders „schone“ Muster ergeben, wiirde ich mich freu- 
en (email: lud.wig@web.de). 
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Anmerkung3: 

Wichtig fur den Eindruck, der beim Betrachten eines Musters entsteht, ist die Grofte der 
„Punkte“, die das 2D-Fenster zeichnet. Entscheidend ist der Fullgrad des „Punktes“, d.h. das 
Verhaltnis der Flache des Punktes zur Einheitsflache des Koordinatensystems. Gute Ergeb- 
nisse erhalt man, wenn der Fullgrad zwischen 25% und 150% liegt. (Mehr als 100% bedeu- 
tet, dass benachbarte Punkte verschmelzen.) Zur schnellen Feinabstimmung kann man das 
2D-Fenster verkleinern Oder vergroftern. 

Anmerkung4: 

Man kann die Funktion g(x) durch ein anderes Polynom mindestens 2. Grades und mit Koef- 
fizienten aus Z[i] ersetzen und weitere Muster studieren. 

Anmerkung5: 

Die Funktion MODCS kann ohne Nachteil durch die Funktion MODC ersetzt werden. Da- 
durch verschiebt sich das Quadrat der Reste urn seine halbe Diagonale. 

Anmerkung6: 

Der Begriff „doppelt-periodisch“ wird tiblicherweise in der Funktionentheorie verwendet. 
Wenn hier eine Funktion „2-dim-periodisch“ genannt wird, soil das bedeuten, dass es zu die- 
ser Funktion zwei iiber Q linear unabhangige Zahlen u und v aus Z[i] gibt, so dass fur alle x 
aus Z[i] gilt: f(x+u) = f(x+v) = f(x). 

Anmerkung7: 

Neben Z[i] gibt es weitere euklidischen Ringe der Art Z[sqrt(k)], namlich zumindest fur k aus 
{3, 2, -1 , -2, -3, -7, -11}. Es ware interessant auch deren Muster zu untersuchen. 


Fragen: 

Fur welche Moduln 

• hat das Muster die Symmetrie eines Quadrats? 

• hat das Muster die Symmetrie eines Rechtecks (bzw. einer Raute)? 

• ist das Muster invariant unter Drehung urn 90°? 

• ist die reelle Oder die imaginare Achse vollstandig besetzt? 

• ist eine Diagonale durch den Ursprung vollstandig besetzt? 

• hat das Muster Unterperioden, d.h. Perioden, die ||m|| / n, n > 2, Punkte umfassen? 


Some comments and questions are following: 

Commentl: 

Interesting patterns can be found for moduli (= divisors) m with: 

a) meZ and odd, e.g. m = 7, 9, 27, 41 , 49, 73. 

b) ||m|| = m m* is a square (i.e. real and imaginary part belong to a Pythagorean triple, e.g. 
m = 40+9i 

Pretty patterns are created by: m = 9+i, 15+i, 39+i, 39+5i, 39+22i, 39+29i, 40+1 1 i, 40+23i, 
40+29i, 40+39i, 41 +i, 41 +9i, 41 +20i, 41 +38i. Patterns for a + b i and b + a i are symmetric. 
For a e Z and a <> 0 (mod 4) seems to be that the patterns with respect to a and a+a-i are 
the same. 

Comment2: 

I would be very grateful to receive any notice about moduli m which give extremely “beautiful" 
patterns (email: lud.wig@web.de). 
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Comment3: 

Important for the impression given by the plot is the point size in the 2D-plot window. I mean 
the ratio of the area filled by one single point with respect to the unit area of the coordinate 
system. Good results can be expected with a ratio between 25% and 150%. (More than 
100% let neighboring points melt into one another.) 

Comments 

Function g(x) can be substituted by another polynomial with order > 2 and coefficients from 
Z[i] in order to investigate or study other patterns. 

Comment5: 

MODCS can be substituted by MODC. The square of the remainders is then shifted by its 
half diagonal. 

Comment6: 

The term “double periodical" is usually used in function theory. We call here a function “2- 
dim- periodical" to indicate that there are two above Q linear independent numbers u, veZ[i] 
such that for all xeZ[i]: f(x+u) = f(x+v) = f(x). 

Comment7: 

Beside Z[i] there are other euclidean rings of kind Z[sqrt(k)], namely at least for k e {3, 2, -1 , 
-2, -3, -7, -11}. It might be interesting to investigate their patterns, too. 


Questions: 

For which moduli m 

• do we receive a pattern with the symmetry of a square? 

• do we receive a pattern showing the symmetry of a rectangle (or a rhombus)? 

• is the pattern invariant under a 90° rotation? 

• is the real or the imaginary axis completely occupied (by points)? 

• is one diagonal passing the origin completely occupied? 

• are there subperiods, i.e. periods which contain ||m|| / n, n > 2, points? 


raWKIOT 
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Histograms of discrete probability distributions 

Benno Grabinger, Guido Herweyers, Josef Bohm 

I received a mail from our friend Benno Grabinger, complaining that it seems to be impossi- 
ble to present the probability function of e.g. a binomial distribution as a histogram. In MS- 
Excel is this no problem at all. 

Benno had an idea to overcome this deficiency by producing a list of many - say 1000 or 
more - outcomes of this distribution and then plot the respective histogram. 

Example: n = 2 and p = 0.5 gives probabilities {0.25, 0.5, 0.25} for x = 0, 1 , 2. Iiste_bin con- 
tains 250 zeros, 500 ones and 250 twos. And this list can be presented as a histogram. 


Lieber Josef, 

mich hat es schon immer geargert, dass 
man von einer Wahrscheinlichkeits- 
funktion die als Liste vorliegt nicht 
ohne weiteres ein Histogramm zeich- 
nen kann (zumindest ist mir das nicht 
gelungen), so wie das z.B. in Excel 
ohne weiteres moglich ist. 


liste_ 1 : =binomPdf(ii,p) 

► { 0.000559,0.005027,0. 

► 


© So kann man die Binomialverteilung nicht als 
Histogramm darstellen. Idee auf der nachsten Seite: 
Erzeuge eine Liste mit den Elementen 0, 1,2, 3, 4, 5, 
wobei jedes Element (ungefahr) so oft vorkommt 
wie der mit 1000 multiplizierte Wert aus der Liste 
binomPdf(5,0.5) 


Jetzt hatte ich eine Idee, wie man das 
doch machen kann. Wenn es dich inte- 
ressiert, so schau dir einmal die Datei 
histo an. Hat diese Idee sonst schon 
jemand gehabt? 

Liebe GriiBe, 

Benno 


Beispiel: fur n = 2 und p = 0.5 ist die 
Liste der Wahrscheinlichkeiten fur x = 
0, 1, 2, natiirlich {0.25, 0.5, 0.25}. 
liste bin erzeugt dann 250 Nullen, 500 
Einser und 250 Zweier. 


0 0.02 0.06 0.10 0.14 0.18 0.22 0.26 0.30 0 . 3 < 
liste_1 



Fur diese Liste kann schlieBlich das 
Histogramm erzeugt werden. 


urliste 

Define m liste| v/_ _xiiO')= 
Func 

Local liste j 


For i, l,diin(>ifl_x) 

liste : =au gment(fi ste, seq (w£_x[i ], K 1, ] ■ 1 0 0 o) ) ) 

End For 
Return liste 
End Func 
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I wrote to our Statistics Guru Guido Herweyers - and he knew the solution: We have to apply 
a “Summary Plot" which is in the German version a „Ergebnisdiagramm“ in order to get a 
histogram of numerical data. Many thanks Guido for your advice. (See also Guido’s contribu- 
tion on Statistics with Tl-Nspire in DNL#85, page 42. I should study the DNLs more accu- 
rately, ©) 

Now I could proceed not only presenting a binomial distribution. I wanted to have all popular 
discrete distributions within one single file working with sliders for the distribution parameters. 

The spreadsheet was filled with the numbers of successes and their respective probabilities. 
For the parameters n (sample size), p (probability of success in geometric and binomial dis- 
tribution), pop (population size in hypergeometric distribution), succ_p (number of successes 
in population - hypergeom. distr.), and finally X (expectation in Poisson distribution) were 
sliders installed. See the spreadsheet and the histogram of the geometric distribution. 
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Benno-Guido-Josef: Histograms of discrete distributions 
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These are the hypergeometric and the Poisson distribution followed by the binomial distribu- 
tion with the approximating normal distribution function superimposed. The sum of all rectan- 
gle areas is always 1. I added some examples which should be solved using the histograms. 




Examples: 

Solve all examples using the graphic representation of the respective clistrubution. 

(1 ) A hunter hits his goal with a probability of 40% for each shot. What is the probability to hit 
his next goal on his fourth attempt? (8.64%) 

(2) Analysis in a rescue organisation finds out that in 47 cases of 1 00 doctors are attending the 
ambulance crafts. Yessterday were 8 missions necessary. What is the probability that three of 
them were attended by a doctor? (24.31 %) 

(3) In a box containing 20 apples are seven apples not of 1st quality. A customer takes 5 of 
them without checking their quality. What is the probability that he buys only 1 st quality 
apples? (8.3%) 

(4) The prob. that a tourist in an certain tropic county will get infected by a tropic desease is 
about 0.0002%. Approximately 30 000 tourists come in per year. What is the probability that 
one of them gets infected? (5.65%) 
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